PweEincThetaAllPhi

 Calculates the incident electric field E from its PWE expression for a fixed r0 but many theta and lambda
 The PWE expression is simply \mathbf{E}=exp(i kM z) \mathbf{e}_x
 The fields Ecr, Ect, Esf given in the results are discussed in the
 supplementary information.

 Parameters:
 - lambda: column vector [L x 1]
           wavelengths in nm
 - epsilonM: scalar or column vector [L x 1]
           epsilon of dielectric of outside region
 - r0:     scalar [1 x 1] 
           spherical coordinate r0 (in nm) of points
 - theta:  possibly row vector [1 x T]
           with spherical coordinate theta of points

 Returns:  stEAllPhi, structure with 3 fields
           containing matrices [L x T]
           representing the three components E_{cr}, E_{ct}, E_{sf} such as: of
           E = E_{cr} cos(phi) e_r + E_{ct} cos(phi) e_theta + E_{sf} sin(phi) e_phi
 - stEAllPhi.Ecr is E_{cr}
 - stEAllPhi.Ect is E_{ct}
 - stEAllPhi.Esf is E_{sf}

 This file is part of the SPlaC v1.0 package (copyright 2008)
 Check the README file for further information

0001 nNbLambda=length(lambda);
0002 
0003 % exp(ikM z) is [L x T], obtained by matrix product of [L x 1] by [1 x T]
0004 phasefact=exp(i*2*pi* (sqrt(epsilonM)./lambda) * (r0*cos(theta) ) );
0005 
0006 % Results are all [L x T] matrices
0007 % They result from Eq. H.16 for e_x and H.76 for E_inc
0008 stEAllPhi.Ecr= phasefact .* repmat(sin(theta),nNbLambda,1);
0009 stEAllPhi.Ect= phasefact .* repmat(cos(theta),nNbLambda,1);
0010 stEAllPhi.Esf= - phasefact;
0011